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Abstract. I review a new treatment of an old idea for light-front quantization of lattice 
gauge theories and give new results from some illustrative calculations: [I] Transverse 
Lattice Gauge Theory; [II] Pure Glue; [III] Heavy Sources and Winding Modes; [IV] 
£>. ■ An Example - Large- N QCD in 2 + 1 Dimensions. 
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<N : I TRANSVERSE LATTICE GAUGE THEORY 

On i 

A Light-Front Co-ordinates 

In the constituent quark model, hadrons are composed of a few quarks moving 
about in empty space and bound by non-relativistic potentials. It is an extremely 
successful model. Before the advent of QCD as the quantum field theory of hadron 
structure, there was little reason to think of hadrons in any other way. But while 
there is consensus that QCD is fundamentally correct, the low-energy dynamics 
of this gauge theory are thought to involve a complicated vacuum with arbitrary 
numbers of relativistic quarks and gluons. How could these pictures be reconciled? 
Indeed, how can one efficiently tackle the problem of relativistic strongly bound 
states at all? 

There is a popular yet poorly understood formulation of quantum field theory — 
Light-Front quantisation — where something like the constituent picture of bound- 
states arises naturally. In fact, it offers much more than that, since it furnishes a 
Schrodinger equation suitable for relativistic many-body systems. 

In the co-ordinate system {x° •j Ob j Ob >j Ob 3 }, where x° is time (I set c = h — 1), one 
usually initializes the wavefunction ^(a; 1 , x 2 , x 3 ) on a hyper-plane a; = constant. If 
the four-momentum of the system is {P°, P 1 , P 2 , P 3 }, the energy-operator P° is the 
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equal-x hamiltonian that evolves the wavefunction in x° according to Shrodinger's 
equation 

This formulation has manifest rotational invariance and is tractable for most non- 
relativistic problems in physics. It is not so suitable when the momenta p t=l > 2 > 3 ~ 
P° (recall that P° contains the rest-energy). 

In light-front (LF) quantisation [1] on the other hand, we (conventionally) define 

, x° ± x 3 , P° ± P 3 
x ± = =- , P ± = ^— . 2 

P ± is canonically conjugate to x T , and we interpret (x~ , P + ) as space and momen- 
tum, while (x + ,P~) as time and energy, x = {x 1 ^ 2 } and P = {P^P 2 } are the 
transverse co-ordinates. The Minkowski metric in the new co-ordinate system is 
such that the invariant length 

ds 2 = 2dx + dx~ - (dx 1 ) 2 - (dx 2 ) 2 (3) 

and x + = x_. The wavefunction ty\ c is initialized on a null hyper-plane x + = 
constant and the LF Schrodinger equation that evolves ^ \ c (x~ , x 1 , x 2 ) in LF time 
x + is 

This equation is now suitable for relativistic problems. For example, \&i c is mani- 
festly Lorentz-boost invariant; a detailed discussion of the advantages and historical 
survey of LF ideas can be found in ref . [2] . 

It is easy to see why LF co-ordinates may help with the constituent question, if 
we consider the energy-momentum relation for a (free) particle of mass \i 

p = ^pr • (5) 

With the conventional interpretation of positive energy (P° > 0) for particles and 
anti-particles, it follows that light-front momentum P + > 0. Because P + is con- 
served and the total P + of the vacuum must be zero (by translation invariance), 
each particle contributing to the vacuum state can have only P + = 0. But accord- 
ing to to (5), these 'zero modes' have infinite light-front energy. If we introduce a 
high energy cut-off, as is always necessary for defining a quantum field theory, these 
zero modes and the vacuum structure they carry will be explicitly removed. Their 
effects on observables below the cut-off should be felt through renormalisation of 
the hamiltonian. If this renormalisation entails quarks getting a constituent mass, 
and the gluon degrees of freedom getting at least a mass gap, if not constituent 



masses themselves, then we have arrived at the basic picture of the constituent 
quark model. 

The trouble with the above idea is that the LF renormalisation procedure is 
still quite poorly understood for theories like QCD, though interesting efforts are 
being made in this direction [8-11]. A formulation of quantum field theory where 
non-perturbative effects normally associated with the vacuum have to appear via 
explicit counter-terms in the hamiltonian is also alien to many people. For these 
reasons there have been some reservations about the consistency of light-front quan- 
tisation. Nevertheless, the ideas are physically appealing and a number of examples 
are known where one can derive the renormalisation resulting from the removal of 
certain non-dynamical zero modes; see Burkardt's lectures at this school in a pre- 
vious year [12]. The attitude I will take in these lectures is that the potential 
applications of light-front quantisation are far too important for it to be simply 
dismissed. I will describe a practical framework for doing calculations in LF QCD 
with a high-energy cut-off, pioneered by Bardeen and Pearson [13,14], which uses 
symmetry to guide renormalisation of the theory. The symmetries that define QCD 
are gauge, Lorentz, and chiral invariance. In these lectures, I will only discuss pure 
gauge theories and infinitely heavy quarks, so chiral invariance will not be treated. 

The method will initially allow all operators in the LF hamiltonian that respect 
symmetries unviolated by the cut-off — allowed operators. After systematically 
truncating to a subset of these allowed operators, one then tunes the remaining 
couplings a posteriori to restore the violated symmetries in observables as best one 
can. This is a physically motivated, systematically improvable framework which 
does not use any fits to experimental data, i.e. it uses only first principles. At 
present, the tests for violations of Lorentz covariance are done in a rather inef- 
ficient way, simply by solving a whole bunch of hamltonians to find the correct 
one. Nevertheless, the recent results obtained in this way [5-7] are surprisingly 
accurate. Once a more efficient treatment of Lorentz covariance is developed, or 
non-perturbative LF hamiltonian renormalisation group ideas are better under- 
stood, many new areas of non-perturbative QCD will be opened up. 

B Transverse Lattices 

Light-front quantisation already has more manifest Lorentz symmetry than any 
other hamiltonian quantisation scheme [1], but typically one cannot avoid breaking 
rotational invariance; the choice x 3 in (2), rather than x 1 or x 2 , is arbitrary. On 
the other hand, we would like to preserve gauge invariance (later this will also 
greatly facilitate the treatment of confinement). This can be done with a lattice 
cut-off [15,16]. In a hamiltonian formulation time remains continous and infinite, 
of course. For a LF hamiltonian P~, whose configuration space is {x~ , x 1 , x 2 }, only 
lattice discretization of the transverse directions x is appropriate. Discretizing the 
longitudinal co-ordinate x~ is not appropriate since it would cut-off large values 
of the conjugate variable P + ; but it's small values of P + that correspond to high 



energy. To remove the high-energy region in a gauge-invariant way we can impose, 
say, anti-periodic boundary conditions on x~ [17]. In the following we use indices, 

v e {0,1,2,3,}, a,f3 e {+, — } and r e {1,2}, and sum over repeated indices. 
We introduce a transverse lattice spacing a and a longitudinal period C. 

To construct an ^i/ (N) transverse lattice (pure) gauge theory on this spacetime, 
we introduce gauge potentials A a in the Lie algebra of SU(N) for the continuum 
directions, and lattice variables for the transverse directions. The particular choice 
of lattice variables which will be convenient for LF quantisation are complex NxN 
matrices M r (x + , x~, x) associated with the transverse link from x to x + ar, at po- 
sition (x + , x~) — colour-dielectric link variables. Because they are linear variables, 
it is easy to identify the independent degrees of freedom. It is sometimes helpful to 
think of them as an average over paths C between these points, with some weight 
p(C), of the short-distance continuum gauge potentials 



Near the continuum limit a = 0, the potentials must change only very slowly over 
many lattice spacings, and M must be forced to lie in the SU (N) group. At larger 
a there is no such restriction however. In fact, for non-abelian gauge theories, it 
makes more physical sense to use disordered complex variables M when the lattice 
cut-off a is quite large. 

Continuum gauge transformations induce the following lattice gauge transforma- 
tions on these variables 



where V G SU(N). The strategy will be to construct LF hamiltonians invariant 
under these gauge transformations and any Lorentz symmetries that have not be 
violated by the cut-offs; this include boosts along a; 3 and discrete Z 4 rotations about 
x 3 . We will seek to enhance the remaining Lorentz symmetries, generically broken 
by the cut-offs, by tuning the couplings of these hamiltonians. 
The continuum QCD action is 



The following gauge-invariant transverse lattice action is the simplest that reduces 
to (9) in the naive continuum limit a — > 0, C — > oo, [14] 
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where 



D a M r {x) = (d a + L4 Q (x)) M r (x) - «M r (x)A a (x + of) , (11) 

A^piaq is the product of link matrices around a transverse plaquette, and U[M] is 
any gauge-invariant potential that forces M into SU (N) as a — > 0, such as 

U[M] = j (Tr{(l-M r t (x)M r (x)) 2 } + (detM-l) 2 ) , (12) 

where A — > as a — > 0. 

More generally, S*[M] can consist of all the combinations of M, F a @, D a M which 
are invariant under gauge and residual Lorentz symmetries. Their couplings are to 
be chosen so as to restore the symmetries violated by the cut-offs. 

II PURE GLUE 

A Colour-Dielectric Expansion 

S[M] contains an infinite number of allowed operators, so to begin calculations 
we need to truncate them to a finite set in some physically reasonable way. Of 
the two cut-offs, a and C, it will be possible to extrapolate the latter. It makes 
sense therefore to restrict allowed operators on dimensional grounds w.r.t. {0, 3} 
co-ordinates. The divergences that appear as C — > oo are of normal-ordering type, 
and thus easily dealt with. We will also demand an interaction-independent (ki- 
netic) momentum operator P + and naive LF parity x + <-> x~ . The main further 
restriction will be to expand the LF hamiltonian P~, that results from 5*[M], in 
powers of M — the 'colour-dielectric expansion'. This only makes sense if, for a 
given a, the P~ that best recovers Lorentz covariance is analytic about M = 
and (classically) minimized there. (This could not be the case in a neighborhood 
of a = 0, where M should be forced into SU(N)). We are not guaranteed to find 
anything at all. There is, however, circumstantial evidence from Euclidean lattice 
work [18] that a Lorentz covariant scaling trajectory in the space of couplings flows 
from the continuum limit to a large a region where M = is the classical minimum. 
Our job is to explictly verify the existence of such a trajectory for the transverse 
lattice by testing this region for signs of Lorentz covariance restoration. 

B Gauge- Fixing 

Light-front quantization is greatly simpled by the LF gauge choice A_ = 0. 2 This 
still leaves ^"-independent gauge transformations unfixed. Suppose we start with 
an action 

2 ) Because we use anti-periodic x~ boundaries, the / dx~ A- zero mode has been removed from 
the theory. Its omission can be shown to renormalise the mass of the M field. 



S = J dx°dx 3 J2 (Tr^D a M r (x) (D a M r (x) 

_-l^Iv{F a/3 F^}-K[M]) , (13) 

where V[M] contains all gauge-invariant products of M up 0(M 4 ) and G 2 — > g 2 /a 2 
as a — > 0. This is still not yet a finite number of operators, but we can obtain a 
finite set by further requiring a degree of transverse locality for operators that are 
products of gauge-invariant operators (see later). 
This action in LF gauge reduces to 

S(A_ = 0) = [ dx°dx 3 Y, (Tr{9+M r (x)9_M r (x) t } + c.c. 

X 

+ ^Tr{(«9_A + ) 2 } + Tr{A + J+(x)} - V X [M]) , (14) 

where 

J+ = i f M r (x) H_ M r f (x) + M r f (x - ar) d- M r (x - of)) . (15) 



We note that the equation of motion for A + is a constraint (no x + -derivatives) 

^- 2 {d-?A + = J + -lTrJ + , (16) 
which can be inverted for A + 

A + = ^(j + -l^J+) . (17) 
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C Canonical LF Quantisation 

The following canonical momenta can be straightforwardly derived in gauge- 
invariant form from the usual definition in terms of the energy momentum tensor 
P a = Jdx-E x T +a ; 

P + = 2 [ dx~Y, Tr{<9_M r (x)<9_M r (x) f } (18) 

X 

P- = J dx- £ V[M] - ^ Trj J+^L J+| + ^ Tr Tr J+ . (19) 

The transverse momenta P and boost-rotations M^ v will be discussed shortly (see 
ref. [2] for a summary of important properties of the Poincare generators in light- 
front formalism). M r is conjugate to <9_M^, and we impose equal-x + commutation 
relations 



M r ,ij{x ,x),(d-M 8tkl (y ,y)) t = -5 u 5 jk 5(x -y )5 x , y <5 ri 



(20) 



where i,j G {1, • • ■ N} are colour indices. 

It is useful to make a Fourier decomposition at x + = with respect to x~ in 
terms of free fields 



M r (x + = 0,x~,x) 
that induces 



l=T% (a- r (k + , x) e-^~ + at (k\ x) e* k+ *~) (21) 
47r Jo vA; + v 



a Ajii (A; + ,x),a* )fci (A; + ,y)] = 5 iJfc 5 jX 5 Xp 5 x , y 5(£: + - k + ) , 
a Ajii (A; + , x), a pM (k + , y)] = , 

A, p G {±1, ±2}, a A ^ = (a A )ji- Acting on a Fock vacuum |0 > defined by 
a x ,ij(k + ,x)\0 >= V k + ,x,i,j,\ , 



(22) 
(23) 

(24) 



°a ij(k + j x ) creates a link-parton on (x, x+aA) on the transverse lattice and carrying 
longitudinal momentum k + . (This mixed momentum-coordinate representation will 
be useful for displaying confinement shortly.) The general Fock space then splits 
into blocks of definite total P + 



{a*(P + )|0 >,a*(P + -k + )a*(k + )\0 >,...} 



(25) 



States of definite P will be constructed shorty. The above Fock basis can be used 
to write down a matrix representation of P~ , which can then be diagonalised to 
obtain LF wavefunctions ^/\ c . 



D Physical States 

The first key point to note at this stage is that, because of positivity k + > 0, 
P~ and P + contain no terms with only creation operators a*. This means that the 
Fock vacuum satisfies : P~ : |0 >=: P + : |0 >= 0. Of course, this is only true so 
long as M is a massive degree of freedom; this is the region of couplings we will 
be exploring for Lorentz covariance. Provided no tachyons appear in the physical 
spectrum of boundstates, |0 > is the physical vacuum of the theory (in the presence 
of our cutoffs). 

Fock space is further simplified by confinement [13]. Since the operator corre- 
sponding to k + is cL, a glance at (19) shows that eliminating A + has introduced 
small- k + singularities. These are cut-off by C, but will blow up as C — > oo for fixed 
a in general. The divergence is avoided if the LF charge Q associated to the current 
J + (15) vanishes 



lim / dx~e ik+x ~J + = 



(26) 



Only certain Fock states satisfy this condition, : Q : |Fock >= 0. They are the 
states invariant under the residual ^'-independent gauge transformations 

M r (x) -> \/(x)M r (x)\/ t (x + of) , d_V = . (27) 

This will include, for example, closed loops of link-partons on the transverse lattice 

Tr{a[(P + -k + ,x)al x (k + ,x)}\0> , 

Tr{4(P + - Jfe+ - k+ - A;+,x)4(A;+,x + aA)al A (A;+,x + 2aA)al A (A;+,x + aA)}|0 > , 
etc.. 

This simple geometrical picture of the finite-energy, gauge-singlets is why we chose 
to work in transverse configuration space but longitudinal momentum space ini- 
tially. 

A remarkable property of the Fock space of connected closed loops of links 3 is 
that it is finite-dimensional in a sector of fixed P + and P for given cut-off C. If we 
write the momentum of the ith parton as 

n-P + 

n = ^, (28) 

where n { G {1/2,3/2, ... , K — 1/2} and K is an integer that plays the role of 
dimensionless cut-off, the Fock space is finite-dimensional because individual mo- 
menta are discrete and bounded and the number of partons is bounded (by 2K). 
This is the basis of the numerical treatment known as DLCQ [19]. Matrices of 
finite dimension can be calculated and diagonalised, and the results extrapolated 
to large K. 



E Momentum Eigenstates 

To test Lorentz covariance it will be important to have explicit forms for bound- 
state wavefunctions at non-zero momentum. So far we have constructed a Fock 
basis diagonal in P + , with modes localised on transverse links. To obtain states 
of definite non-zero transverse momentum P, we can take a gauge-singlet p-link 
shape, say, and translate it over the transverse lattice sites y according to 

£ e* p ^ + *) Tr{a Ai Xl + y ) ^ x 2 + y) • • • a{ p {k^ x p + y)} |0> . (29) 
y 

We must have 



3 ) These are the relevant ones in the N = oo limit for example. 



1=1 

Xj = Xj_! + a Aj_i , i<l<p (30) 

i=l 

and it is convenient to set the phase convention by adopting a momentum-weighted 
'centre-of-mass' 

*=^XX(** + ^) . (31) 

This is consistent with the canonical gauge-invariant form of the M_ r component 
of the boost-rotation tensor 

M_ r = 2 f dx~ J2(%r + ad rs /2) Tr{9_M s (x)9_M s (x) t } (32) 

X 

which gives 

Q -WM. T fl t x) e i6"M_ r = a f ( £+ x) e -ifc+b.(x + aA/2) (33) 

e -i^M_ r |^ c(p+) p)) = |^ ic(p + p _ bp+)) ? (34) 

Thus, non-zero P states in the Brillouin zone (— 7r/a,7r/a) can be conveniently 
obtained from those at P = 0. 4 Longitudinal boosts generated by 

M_+ = 2 f dx-J2%~ Tr{9_M r (x)9_M r (x) t } (35) 

X 

are even more trivial; up to unimportant global phases they simply rescale longi- 
tudinal momentum 

e isM -+ |M> lc (P + , P)> = \^ c (e~ s P + , P)) . (36) 

Thus, longitudinal momentum fractions Xi = kf / 'P + are Lorentz-boost invariant. 
A perfectly relativistic state should have a dispersion of the form 

M 2 = 2P + P- - P 2 (37) 

and so diagonalising P~ in a basis of fixed P + and P is equivalent to finding the 
masses M. of boundstates. Because of the transverse lattice however, there will 
be corrections to (37) (see later) which we will try to minimize as part of the 
optimization of Lorentz covariance. 

4 ) The construction of the operator P itself is more problematic however, since the canonical 
expression is not gauge- invariant, while the simplest gauge-invariant extensions are not quadratic. 



Ill HEAVY SOURCES AND WINDING MODES 



Pure QCD is characterised by a single scale in terms of which all dimensionful 
quantities, such as hadron masses, can be expressed. A convenient measure of this 
scale is the string tension a, the asymptotic slope of the linearly-rising confining 
potential between heavy sources. Moreover, by measuring this slope in both trans- 
verse lattice directions x and the continuum space direction x 3 , and requiring it to 
be rotationally invariant, we fix the lattice spacing a in dimensionful units. 

A Winding Modes 

In the transverse directions, a more accurate measurement of a can be made by 
examining winding modes rather than heavy sources. By making the transverse 
lattice compact in direction A say 

x = x + a\D\ , (38) 

where D\ is the number of transverse links in direction A, we can construct a basis 
of Fock states that wind around these directions, e.g. 

Trja^ x) 4 2 (k+, x + aA x ) • • • a[ p x + a\D x - a\ p ) } |0> (39) 

has winding number 1 in direction A. The mass spectrum of such winding modes 
should rise linearly with D\ with slope a as D\ — > oo. Unlike the potential between 
heavy sources, there are no 'endpoint' effects, and so the asymptotic linear rise 
should set in more quickly, especially for the lowest mass eigenvalue 

M — »■ o- T |n|a , |n| — > oo , (40) 

where n = (Di^i, D 2 n 2 ) and rii, n 2 are winding numbers in directions (1,0) and 
(0, 1). (Note that we cannot construct winding modes around a compactified x 3 as 
this clashes with the choice of LF coordinates). The suffix T indicates a transverse 
measurement of a. 

B Heavy Sources 

To measure the potential between two heavy sources in transverse lattice gauge 
theory [20], we will start with a heavy scalar field 4>(x + ,x~ of large mass p in 
the fundamental representation of SU(N). In addition to the pure glue 'link-link' 
interactions, we must now also consider interactions with 0's. The simplest are 



{d a <fi - i^A a ) (d a <j) + iA a <j>) , -p 2 0V 



(41) 



such that <fi couples to M via A a . 5 In LF gauge A_ = this leads to the substitution 

J + - ^p + urc glue [M] + JLrcc [0] (42) 

(JLrceh = -id-M*j + ifad-tf (43) 

Now for some kinematics: Let P f ° u represent the full 2-momentum of a system 
containing h heavy particles. It is convenient to split the full momentum into a 
"heavy" part plus a "residual" part P a due to interactions, 

PU = phv a + P a , (44) 

where v a is the covariant velocity of the heavy sources, v a v a = 1. Note that P a , as 
defined here, is not positive definite. The full invariant mass-squared (at P = 0) is 

M 2 = 2P+ u P full = (hp) 2 + 2hpv + p- + 2P+ (P~ + hpv~) (45) 

The choice of v + is arbitrary and it is convenient to choose it such that P + = 
0. Consequently, v + P~ is just the shift of the full invariant mass M. due to the 
interactions: 

M = hp + v + P- + O (l/p) . (46) 
Thus, v + P~ is the usual energy associated with the heavy quark potential. 



C Fock Space 

We define creation-annihilation operators associated with the heavy field: 

U* + = 0, *- x) = /I - / 0== ( H k, x) e-^>- + d*(k, x) — " ) 

(47) 



bi(k,x),bXk,y) 



: , y 5(k — k) , etc. 



(48) 



The e tpVaXa term removes an overall pv a from the 2-momentum. b*(k,x) creates a 
source particle at site x carrying residual momentum k, while d* creates an anti- 
particle. The residual gauge invariance in A_ = gauge again leads to confinement 
into singlet states. An example of a gauge-invariant Fock state with two co-moving 
sources at x and y respectively, maintaining a fixed a x 3 -separation L is 



f 

j — i 



die 



6*a-Piink/2,x) a^x 



t (£; p ,yKH-P link /2,y)|0> (49) 



5 ) The full set of operators allowed at leading order of the colour-dielectric expansion is described 
in rcf. [7], but we skip the details here. 



where X)f=i — -Pimk- To obtain LF eigenfunctions of P~ we must take linear 
combinations of all possible Pn n k > 0. In practice this requires us to introduce a 
high-energy cut-off on Pn n k, in addition to the DLCQ one, K. v + P~ is then a finite 
dimensional matrix (at least for the connected Fock states that dominate in the 
large- iV limit) whose eigenvalues may be extrapolated to infinite cutoff. 

The lowest eigenvalue of v + P~ can be compared with the popular phenomeno- 
logical form of the potential between heavy sources of spatial separation R 

V(R) =aR + Cl + ^ (50) 
R 

R = \Ja 2 \n\ 2 + L? Source Separation, 

where an = y — x and L = y 3 — x 3 . For generic couplings in P~ the eigenvalues will 
not simply be a function of R, i.e. a, c±,C2 extracted from (50) will depend upon 
L and \n\ separately. Once rotational invariance is restored, the value of a can be 
set from experiment by the masses and decays of heavy mesons using V(R) in a 
conventional Schrodinger equation. A typical value is yfo ~ 440MeV, though one 
must remember that this will includes effects of light-quarks which are neglected in 
pure gauge theory. Higher eigenvalues of v + P~ should correspond to hybrid heavy 
mesons. 



D Setting the Scales 

We now have everything necessary for performing a first-principles calculation. 
Let us now combine results from each sector — pure glue, winding, heavy sources. 
Pure Glue 

For generic couplings in P~ , the glueball eigenstates at fixed momenta (P + ,P) 
have eigenvalues which we may expand in powers of transverse momentum thus 

2P+P- = G 2 N(M 2 + M 2 1 a 2 \P\ 2 + 2M 2 1 a 2 P 1 P 2 + Mla 4 \P\ 4 + ---) . (51) 

Here we have chosen to factor out the coupling G 2 as an overall mass scale (G has 
units of energy). Ai Q ,/Ai,M 2 , • • • are then dimensionless numbers which we cal- 
culate when diagonalising P~ . Eq.(51) is not, of course, in the form of a relativistic 
dispersion relation (37) in general. The coefficients, however, are functions of the 
couplings at our disposal and may be adjusted to regain a relativistic form. 
Winding Modes 

Let us write the lowest winding eigenvalue (40) as 

2P + P- = a 2 a 2 \n\ 2 = G 2 iV>V|n| 2 , |n| -> oo . (52) 

Again, we factor out G 2 to set the scale, and W is a dimensionless measured 
quantity dependent on the couplings. 
Heavy Sources 



Similarly, at fixed sources separation R 



v +p- = a R = G 2 NSR , R -> oo . (53) 

S is the dimensionless measured quantity that we can adjust with the couplings in 
P~. 

The following conditions are necessary for Lorentz covariance of eigenstates 

a T — a — constant (54) 
Ml a 2 G 2 N - 1 = (55) 
ATi = (56) 

Of course, there are further conditions that one could impose on solutions. But as 
we shall see in the next lecture, in practice these seem to be the most significant 
when trying to remove lattice discretization errors to a first approximation. 

To test for (55) we need the dimensionless combination a 2 G 2 N = W/S 2 that 
follows from (52) (53). Combining (51) (53) allows us to express boundstate masses 
in units of the string tension 

2P+P-(P = 0) = M 2 = ^f- (57) 

and also the lattice spacing itself a = JW/Sa. 



IV EXAMPLE - LARGE- TV QCD IN 2 + 1 DIMENSIONS 

We need to do a test calculation that can be accurately compared with known 
results and is relevant to non-abelian gauge theories in 3+1 dimensions. SU (oo) 
gauge theory in 2 + 1 dimensions provides an ideal test. It is physically very similar 
to SU(3) pure gauge theory in 3 + 1 dimensions, having dynamically-generated 
linear confinement and a discrete spectrum of glueballs; it is even quantitatively 
similar. But factorization of colour Traces in the Fock space states in the large- 
N limit and having only one transverse dimension make an accurate calculation 
feasible. Moreover, good Euclidean lattice Monte Carlo (ELMC) data has recently 
become available for this problem [21], which we can take as 'the right answer'. 6 
The results presented here have been obtained with B. van de Sande and are new 
or improvements on our previously published work [5]. 

The discussion of the previous lectures is trivially adapted to 2 + 1 dimensions. 
Here x ± = (x° ± x 2 )j\pl etc. and x 1 is the single transverse lattice dimension. 
Up order order M 4 , we have in the large N limit with G 2 N finite (suppressing 
transverse indices for clarity) 

6 ) The calculations of Teper are based on the traditional Wilson Euclidean lattice path integral 
formulation of gauge theory. The large N calculations were in fact performed about the same 
time as the light-front work, and to some extent motivated by it. 



p- = I dx- - x Tr { J+ ^ J+ } + w TrJ+ w TlJ+ + ^ 2Tr { MMt } 

+ — Tr{MM t MM t } + — Tr{MMM*M*} + (Tr{MM f }) 2 (58) 



where, in anticipation of always taking the C — > oo limit, we have also used power- 
counting in longitudinal coordinates to limit the number of terms. There are in 
principle an infinite number of different products of individually gauge-invariant 
operators (like the A3 term), obtained by separating the operators in the transverse 
direction. We have also assumed a degree of transverse locality by keeping only the 
most local product to a first approximation. This approximation, like the colour- 
dielectric expansion in gauge-invariant powers of M, can of course be systematically 
investigated. 

To maintain a trivial LF vacuum, we shall be considering the region ji 2 > with 
tachyon-free physical spectrum. G has dimensions of energy and will be used to 
express dimensionful quantities; it is convenient to form dimensionless versions of 
all the other couplings 

m2 = ^v' k = ^k- (59) 

The basic technique we follow is to search the space {m, h, I2, h} for a one param- 
eter trajectory on which observables show enhancement of Lorentz covariance — a 
'Lorentz trajectory'. Moving along this trajectory should correspond to changing 
the spacing a, eventually taking us to the continuum limit. We will, in fact, be 
prevented from reaching a = by our self-imposed restriction m 2 > 0. Neverthe- 
less, assuming gauge and Lorentz covariance define pure QCD, observables should 
be invariant along the (exact) Lorentz trajectory. 7 Of course, we can only obtain 
some approximation to such a trajectory with a finite number of couplings. 

To get some intuition for the space of couplings, let us construct charts that 
display the 'size' of violations of Lorentz covariance. We can use a x 2 -test based 
on the variables (54))(55) (condition (56) is absent in 2 + 1-dimensions) . Condition 
(55) in fact provides a number of variables, one for each boundstate (glueball) in 
the low-lying spectrum. To begin with, let us set the lowest glueball mass to the 
known value M. = 4.05-^/0" [21], and try to get an isotropic speed of light (55) 
for this and the heavier low-lying glueballs. With a particular choice of variance 
assignments in the x 2 -test, we can make plots of {x^ in , m 2 , /«} for each /j, where 
Xmm 1S ^e minimum value of x 2 with respect to the couplings except m 2 and /«. 
For technical reasons 1% almost completely decouples (we can set it to almost any 
fixed large value, see ref. [5]) so we only display charts for Z x (fig.l) and l 2 (fig- 2). 

They each show an unique, narrow valley, running from small to large m 2 , at 
the bottom of which Lorentz covariance is optimised. One finds that changing the 



7 ) The idea is similar to the renormalised trajectory in the renormalisation group associated 
to a fixed point with one relevant direction. However we are not performing any explicit RG 
transformations here! 
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FIGURE 2. x 2 -chart for l^. Darker shades correspond to lower x^i 
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FIGURE 3. Variation of the transverse lattice spacing along the Lorentz trajectory. The fit is 
1.275m + 1.23. 

details of the x 2 test changes these charts somewhat, except in the neighborhood of 
the valley in each case. Thus one finds a degree of universality in the results. If there 
were no Lorentz trajectory, one would not expect to obtain a robust valley. Having 
roughly located the candidate Lorentz trajectory, one can increase the resolution 
in its neighborhood, and perform more efficient iterative searches for the bottom 
of the valley. 

In a first-principles calculation we cannot use the result for the lightest glueball 
mass in o units — this should be predicted by the method. Performing a new 
X 2 search in the neighborhood of the valley, but now including a variable to test 
(54), one finds a similar position for the valley bottom, without having to input 
the lightest glueball mass. The x 2 P er degree of freedom is always about 1 near 
the Lorentz trajectory, as it should be for reasonable criteria. Results as one moves 
along the valley bottom are displayed in fig. 3 and fig.4. We see that the lattice 
spacing gradually decreases with m 2 as one moves along the Lorentz trajectory, 
but never becomes zero for m 2 > 0. Despite fluctuations, the lightest glueball mass 
scales in a way roughly consistent with the correct continuum answer. 

Although the lightest glueball (a J PC = ++ ) is behaving covariantly along 
the Lorentz trajectory — the speed of light deduced from the ++ 's dispersion is 
isotropic to with % 2-3 — fluctuations are still present in fig.4. These are mostly due 
to the difficulty in accurately establishing the scale a. Most of the fluctuations tend 
to cancel if we consider ratios of glueball masses, and good scaling is obtained along 
the Lorentz trajectory [5]. For higher glueballs with nearly-covariant wavef unctions, 
such as the and its excited state 0~~, this provides rather accurate mass ratio 
determinations in the large- N limit; see figure 5. Finite- N results from conventional 
ELMC in ref. [21] were fit to A + B/N 2 , at confidence levels of order %50 - 70. 
Including our N = oo data, and fitting to A + B/N 2 + C/iV 4 , improves this to %95 
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FIGURE 5. Variation (or lack of it!) of the glueball mass ratios /0++ and 0~/0++ with 
N. The fits are 1.349 + 2.9U/N 2 - 13.968/iV 4 for 0~ and 1.824+ 1.274/7V 2 - 7.US/N 4 for 0"". 
The error on the N = oo results is from extrapolation to C = oo. 
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FIGURE 6. The heavy-source potential. Solid line is fit to potential for sources with 
x 2 -separation only; data points are values at one-link transverse separation and x 2 -separation 
Ly/C^N = 0, 2.5, 5. 

for 0~/0 ++ and > %99 for (T~/0 ++ glueball mass ratios! 

The heavy-source potential calculated at a point on the Lorentz trajectory is 
displayed in fig. 6. It shows good restoration of spatial symmetry. The potential in 
the continuum spatial direction x 2 is a fit to 

o 1 78 

v+P~ = 0.15ALG 2 N + 0.183VG 2 N - —— (60) 

Lj 

at the point on the Lorentz trajectory of lowest overall \ 2 ■ One must be careful 
when interpreting (60) since the Coulomb potential in 2 + 1 dimensions is logarith- 
mic. The form (60) should be appropriate except at the very smallest L, where 
Coulomb corrections are expected. The 1/L term is a universal correction expected 
on the grounds of models of flux-string oscillations [22]. Universality implies that 
its coefficient should be invariant along the Lorentz trajectory. In reality, we find 
that it drifts slowly, a symptom that our approximation to the Lorentz trajectory 
is not an exact scaling trajectory and/or the form (60) is not sufficient to fit the 
potential. The coefficient 0.178 at the lowest x 2 is nevertheless close to the theo- 
retical value dn/24, where d is the number of transverse dimensions in which an 
ideal thin-string can oscillate (d — 1 here). But it is larger by 40% (0.178 gives 
d ~ 1.4). We found the same excess by another method to measure the 'central 
charge' d + 2, using the density of glueball states [4]. In ref. [4] we observed that 
this excess in the central charge comes from longitudinal degrees of freedom in the 
large- N QCD flux string, a mode of oscillation which an ideal thin-string does not 
possess. 
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FIGURE 7. Variation of glueball mass ratios to the lightest mass (0 ++ ), as more operators are 
included in the hamiltonian. |c — l| av measures the average deviation from 1 of the speed of light 
in the transverse direction, for the displayed glueballs. 

One may wonder whether results really do improve as more operators are added 
to the hamiltonian. For a given \ 2 test, adding more operators to the hamiltonian 
will inevitably bring one closer to the Lorentz trajectory as measured by this test. 
We see from fig. 7 that the masses of glueballs, which are not part of the test, 
rapidly approach the correct values as more operators are added and covariance is 
improved. 8 



A Future Work 

The transverse lattice idea is already quite old [13]. The breakthrough, that has 
been made over the last year or so, was to show that in 2 + 1 and 3 + 1 dimensions 
large- N gauge theory exhibit a unique Lorentz covariant scaling trajectory on coarse 
transverse lattices [5-7]. Glueball masses on this trajectory are consistent with 
known results. Thus there is good reason to believe that the Lorentz-invariant 
wavefunctions obtained with this method are accurate also. They are completely 

8 ) The 2 h state appears to overshoot the correct mass. It is the only state whose wavefunction 
actually becomes less covariant when Ai is turned on. 



new results, essentially unobtainable within any other quantisation scheme. It will 
be interesting to see what implications they have for experiment. The next step 
towards this goal is to couple these pure gauge theories to propagating quarks 
[13,23], and then enforce the Lorentz and chiral symmetries that define QCD to 
determine the couplings of the hamiltonian. We also need to understand better how 
to extrapolate results into the small a region. This is important not only formally 
for the existence of the continuum limit, but to compare calculations of hadronic 
structure with hard process experiments where the usual factorization theorems 
simplify the analysis. We are only at the beginning of the development of this 
subject, and it is difficult to know how far it will progress, but there does appear 
to be too much truth in it to ignore. 
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